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The prototypical phase change material GeTe shows an enigmatic phase transition at Tc ~650 
K from rhombohedral {R3m) to cubic (FmSm) symmetry. While local probes see little change in 
bonding, in contrast, average structure probes imply a displacive transition. Here we use high energy 
X-ray scattering to develop a model consistent with both the local and average structure pictures. 

We detect a correlation length for domains of the R3m structure which shows power law decay 
upon heating. Unlike a classical soft mode, it saturates at ~20 A above Tc. These nanoclusters 
are too small to be observed by standard diffraction techniques, yet contain the same local motif 
as the room temperature structure, explaining previous discrepancies. Finally, a careful analysis 
of the pair distribution functions implies that the 0.6 % negative thermal expansion (NTE) at 
the R3m -Fm3m transition is associated with the loss of coherence between these domains. 

PACS numbers: 64.60.Cn, 61.05.C 


Introduction: Materials which can be rapidly switched 
between states with contrasting properties are useful for 
memory applications [1]. Examples include the chalco- 
genide phase change materials which are now making 
their way into consumer devices. These materials can 
be switched between crystalline and amorphous states 
using laser heating on a nanosecond timescale [5]. The 
large change in conductivity, combined with kinetic 
trapping of the amorphous state, allows them to be 
used as two-state bits. Understanding the mechanisms 
of phase transitions in these materials is thus of key 
importance. However, a cursory review makes it clear 
that their simple structures and compositions belie 
structural complexity. When doped with antimony, 
GeTe produces one of the most promising phase change 
materials, yet this simple parent phase is still not fully 
understood despite decades of research. In common 
with many other pnictides and chalcogenides, strong 
chemical constraints on bond lengths are found. It is 
this interplay between covalent bonding, coherent lattice 
vibrations and electronic properties which makes these 
materials so interesting. In particular, GeTe shows a sec¬ 
ond order Jahn-Teller (SOJT) distortion j3] of the first 
coordination sphere, due to the formation of Tellurium 
lone pairs. At room temperature, this is manifested as 
an ordered arrangement (Fig. la) of three long and three 
short Ge-Te bonds in each GeTee octahedra, this causes 
a rhombohedral distortion. Upon heating however, the 
rhombohedral angle decreases in an order parameter-like 
manner [4], and a structural phase transition is found 
near 650 K. The resulting cubic rock salt structure 
imposes a single symmetry equivalent Ge-Te distance. 
Average structure probes thus imply that the transition 
is of second order. Spectroscopic measurements lend 


weight to this interpretation, for example inelastic 
neutron scattering and phonon calculations show 
that a triply degenerate phonon condenses at Tc. 
Similar observations have been made with Raman 
scattering [^. In great contrast however, pioneering 
local structure measurements [7] using EXAFS showed 
that the Ge-Te bond lengths are essentially unchanged 
upon heating up close to the melting point. In fact, the 
local SOJT distortion even persists into the liquid phase 
[5]. Unifying the local and averages structure in GeTe 
is thus an important and thus far unsolved problem. 
Here we use high-energy x-ray scattering to address 
this question through so-called total scattering. This 
technique utilises both the coherent Bragg diffraction 
as well as diffuse scattered x-rays, to calculate a pair 
distribution function (PDF) and is thus sensitive to 
both local and long range order. This technique is 
also more sensitive than EXAFS to intermediate range 
order on the nanometer length scale. One earlier 
report applying the PDF technique to GeTe exists [3]. 
This short work confirmed the presence of two Ge-Te 
distances in the high temperature cubic phase. In this 
Rapid Communication, we re-visit GeTe as a function of 
temperature, and apply full structure refinements over 
a range of length scales. We show that: i) the 650 K 
transition corresponds to the loss of coherence between 
nanoscopic distorted domains. These are too small (^ 
20 A) to be observed by diffraction methods, and persist 
close to the melting point, ii) We show that the negative 
thermal expansion |10j at the transition is likely driven 
by increased dynamical disorder, as the average atomic 
distances distances within the distorted domains are 
essentially unchanged. Experiment : Our sample of 
GeTe was the same as used in the original neutron 
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FIG. 1: (color online) a) Structure of R3m GeTe. The pur¬ 
ple spheres represent the Ge atoms, while the golden spheres 
represent Te atoms. One GeTeg octahedra is shown in cut¬ 
away, emphasising the off-centre distortion which generates 
a long and a short Ge-Te distance; b) Refined reduced unit 
cell volume for GeTe as a function of temperature. The open 
symbols are results from our Rietveld refinements, while the 
closed symbols are from real-space refinements against exper¬ 
imental pair distribution functions in the range 30< r < 50 
A as described in the text. 


diffraction investigation [4]. A small piece of the crystal 
boule was ground into a fine powder for the measure¬ 
ments reported here. The sample was sealed in a quartz 
capillary and diffraction experiments were performed 
on the ID15B high energy x-ray scattering beam line of 
the European Synchrotron Radiation Facility (ESRF), 
in Grenoble, France. The incident energy was 87 keV, 
and the scattered x-rays were detected with a Mar345 
image plate. This was placed close to the sample to 
obtain data reaching high momentum transfers (~30 
A“^). At each temperature point, multiple frames were 
collected and averaged in the close configuration, before 
the detector was translated further away to collect a 
data set suitable for conventional Rietveld analysis as 
in previous investigations m- Two dimensional data 
were azimuthally integrated [12] using pyFAI, and the 
PDFGetx3 and PDFgui packages were used to ex¬ 
tract m and model El the pair distribution functions 
respectively. The structure factors were truncated at 
Qmax=‘26 before Fourier transformation. Rietveld 



FIG. 2: (color online) (a) Nearest-neighbour Ge-Te bond dis¬ 
tances extracted from Rietveld refinements (open circles) and 
Gaussian fits (solid circles) to the pair distribution function 
data; b), c). Nearest-neighbour peak splitting observed in the 
pair distribution functions in the R3m and Fm3m phases. 
The lines are the results of fitting two Gaussians and a slop¬ 
ing background. 


analysis was performed El US] using GSAS and the 
EXPGUI interface. 

Results : Turning first to the results of the Rietveld 
analysis, in agreement with earlier reports, a clear rhom- 
bohedral splitting was observed at room temperature. 
Upon heating to ~ 650 K, the splitting reduced until 
a metrically cubic lattice was obtained m HU] . Also 
in agreement with earlier results, the reduced unit cell 
volume shows a ^0.6 % collapse at the R-G transition 
(Fig. lb). As a further check on the high temperature 
symmetry, we also collected a data set in the cubic phase 
using the crystal analyser diffractometer ID31, also at 
the ESRF. This instrument m has an approximately 
100 X better angular resolution than our in-situ data. 
The c-GeTe peaks were essentially resolution limited 
(not shown). This rules out lattice microstrain in 
the cubic phase, and indicates a dynamic origin for 
the nanoscale inhomogenieties inferred from the data 
presented below. The refined values of the symmetry 
inequivalent Ge-Te bond lengths are shown in Fig. 2a 
as a function of temperature. These can be seen to 
converge on a (symmetry imposed) average value at Tc, 
in a manner consistent with earlier work. 

The PDFs which we have calculated from our high 
energy x-ray scattering data explicitly contain the diffuse 
scattering signal, not just the information encoded in the 
Bragg reflections. Upon Fourier transforming our data, 
it is immediately apparent that a long and a short Ge-Te 
distance persist to the highest measured temperatures 
(Fig. 2). A simple model-independent fitting of the 
first coordination sphere bond lengths from the pair 
distribution function data highlights the dichotomy 
between local and average structure measurements. As 
shown in Fig. 2b and 2c, two contributions to the Ge-Te 
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distances can be distinguished at all temperatures. 
Gaussian fits as a function of temperature show that 
these only slowly converge, with an estimated transition 
temperature well above the melting point (~1700 K). 
The resultant energy scale (0.15 eV), is consistent with 
strong covalent bonding m, and explains why the 
SOJT distortion persists even into the molten state [5]. 
Clearly this degree of freedom is of secondary relevance 
to the mechanism of the ferroelectric transition, which 
has a characteristic energy scale of around 55 meV. 

The above results contain little information beyond 
that contained in previously published works. In order 
to go further, we must establish the sensitivity of our 
local structure measurements to the average structure 
as measured by e.g. diffraction. It is well established 
that structure refinement against the pair distribution 
function should tend to the average structure at larger 
values of r in real space m [18]. This is due to the 
Fourier relationship between reciprocal and real space. 
Crudely put, we can reduce the sensitivity to the broad 
diffuse scattering by choosing an appropriate frequency 
window, delimited by low-r and high-r cutoffs. We 
therefore performed real space refinements of the GeTe 
crystal structure in the range 30 < r <50 A as a function 
of temperature. As described elsewhere m, we used a 
damping parameter calculated using a Ce 02 standard to 
account for the intrinsic instrumental resolution. This 
yielded the unit cell volumes shown using solid symbols 
in Fig. lb. This nicely reproduces the Rietveld result, 
proving that the data content of our PDFs extrapolates 
between the local and average structure limits. 

To recap, the PDF of GeTe: i) Shows essentially no 
change in nearest-neighbour distances at the structural 
phase transition and ii) Tends to the known average 
structure result at longer distances. The key question 
is therefore what happens on intermediate length 
scales upon warming through the ferroelectric phase 
transition. To answer this question, we performed r- 
dependent refinements of the R3m crystal structure for 
all temperatures. This model was found to describe all 
data sets satisfactorily at low-r. By plotting the Ge-Te 
bond-lengths as a function of r-max and fitting the 
spherical domain envelope function [I8j . we extracted 
a correlation length for the local order. Note that the 
fitting was done in a self-consistent manner, using a 
single variable to describe the correlation length for 
both bond lengths. Close to room temperature, the 
correlation length is essentially infinite, yet at higher 
temperatures, it begins to decrease. As shown in Fig. 
3c, the correlation length, decays as ^ ex . This is 
reminiscent of the soft-modes observed experimentally 
with neutron and Raman [5] spectroscopies, which 
also show power law behaviour near the transition. Our 
observations are strikingly different however, as the 
correlation length does not tend to zero, but reaches 
a constant level of ^ ~20 A in the cubic phase. This 



FIG. 3: (color online) (a) Real-space dependence of the near¬ 
est neighbour bond distance as htted to the pair distribution 
function at 598 K. Lines show fits of a spherical particle en¬ 
velope used to extract the correlation length; b) As in panel 
a), except data shown are for the Fm3m phase at 798 K; c) 
Extracted correlation length as a function of temperature. A 
power law fit is shown in the R3m region. The inset shows 
the power law scaling as a function of reduced temperature, 
r = , where is 650 K. 


is essentially unchanged until the highest temperature 
we have measured (900 K), and we speculate that 
these nanoscale domains persists all the way up to the 
melting point of GeTe at 998 K. The melt is known to 
contain |H] the same local SOJT distortion as found at 
all temperatures here, and thus may well correspond to 
the loss of coherence in the ~20 A domains discovered 
herein. Note that a simple analysis of the power law 
scaling (inset Fig. 3c) provides a critical exponent of 
v' ^1.7(1). This is far from the mean field value of 
and shows that the transition is far from a conventional 
soft-mode driven transition. 

Finally, we wish to comment on the negative thermal 
expansion observed at the structural phase transition 
(Fig.l b). This phenomena can be driven by many 
diverse mechanisms, such as gross changes [20] in 
electronic structure, or purely by low-energy vibrational 
degrees of freedom. The latter is common in framework 
materials with low coordination numbers |2IL 122] . yet 
is more difficult to justify in a close-packed mate¬ 
rial like GeTe, which also shows abrupt, rather than 
continuous NTE. While our results do not provide a 
definite mechanism for the volume collapse, we can 
associate it with the 2 nm distorted domains discovered 
above. A colour map of the low-r region of the pair 
distribution functions is shown in Fig. 4a as a function 
of temperature. No major anomalies can be observed. 
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FIG. 4: (color online) a) Temperature dependence of the low- 
r pair distribntion function of GeTe, several important dis¬ 
tances are highlighted. The dashed line shows the tempera¬ 
ture at which the bulk R3m to Fm3m transition occurs; b) 
Temperature dependence of peak positions in the pair distri¬ 
bution function. Note the crossover to behavionr mirroring 
that of the nnit cell volnme at higher values of r. Lines are 
guides to the eye, and error bars are shown, bnt are often 
smaller than the symbol sizes. 


which might indicate a reduction in atomic distances. 
In fact, model independent fitting (Fig. 4b) shows that 
e.g. the Ge-Ge/Te-Te peak at 4.1 A smoothly moves to 
high-r. We performed this analysis at several positions, 
attempting to find peaks with low multiplicities. This 
becomes challenging as the length scale increases, as 
the number of distances probed in the PDF is oc 47rr^. 
Nevertheless, clear signatures of the volume collapse 
are apparent as the length scale increases to ^ 20 A, 
matching the size of the distorted domains discovered 
above. We note that tellurium-based alloys generally 
show NTE even in the liquid state. This has been 
convincingly attributed to the interplay between entropy 
and volume |23j . Although our measurements are a time 
average of the thermally occupied states in the system, 
and we are a priori unable to distinguish between static 
and dynamic disorder, we speculate that the latter is 
equally important in ’crystalline’ GeTe. 


In summary, the above results show that the structural 
phase transition found at 650 K in the prototypical phase 
change material GeTe is of an exotic nature. The energy 
scale relevant to the nearest-neighbour SOJT distortion 
is apparently relevant only at temperatures far exceed¬ 
ing the melting point. In contrast, a coarse-grained 
order-disorder model, made up of 2 nm domains, can 
explain the results of both local and average structure 
probes. Future measurements sensitive to the length 
scales probed here, as well as dynamics, will be use¬ 
ful to further understand this important parent material. 
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